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SUMMARY 


The working fluid in a turbcmachine stage can excite vibration of the 
^oto^ must be included as an energy source for the system. The system 
then corprises both fluid and rotor structure, and in stating a stability 
criterion the daitping of both must be considered: 6 =6 +6 

> 0. total fluid rtech 

Mechanical vibration perturbs the passage flow conditions v\hich then 
excites or danpens the inducing vibration. Positive excitation (negative 
daitping) results from tee positive cross-coi:pled coitponents of an unbalance 
t>y the vibratory displacement as well as the orbital velocity 

of the rotor. 

The equation of motion of flow in a turbomachine shews that unsteady 
flew is necessary for work transformation. The anplitude, frequency and 
ptese angle of the unsteady flow phenomena must then be found to coiipare 
with the structural dynamics of the turbomachine. 


Unsteady flew is identified by frequency and size of participating 
fluid ensemble. There occurs at blade passing frequency axial eorpressor 
and turbine blade flutter , and response to blade wakes and condensation 
shock wave oscillations. Near shaft rotational speed inlet distortion, 
gusts, rotating stall as well as self-excited shaft displacement and orbital 
velocity dependent forces are imposed upon the rotor. Surge occurs at 
system resonance frequency. 


The quasi-stea(^ conputer analysis of the perturbed centrifugal im- 
peller ^ssage flew is reviewed. The fluid danping coefficient 6 . , 

linear in the orbital velocity, is defined and based on 115 stage^^^^'^ 
calciiLations, the average total daitping coefficient per stage needed for 
stability is 6 ^ t nr- 


total 


> 1.85. 


* Currently with Borg-Wamer Research Center, Des Plaines, Illinois. 
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THE SYSTEM STABILITY CRITERION 


There is a growing awareness of the fact that the flow of the working 
fluid in a turbcmachine stage can interact with the structure of the turbo- 
machine rotor in such a way that the distribution of the working fluid 
pressure over the passage surfaces is perturbed by the lateral vibration of 
the rotor due, for exaitple, to residual mass unbalance. As a function of 
both flow conditions and passage gectnetry, the pressure perturbaHons , when 
sunroed over the surfaces, itay result in unbalanced forces \^ich introduce 
energy into the rotor and enhance the vibration. Such a vibration is 
deemed self-excited. On the other hand, the flow and geometry conditions 
may result in a net reduction of energy in the rotor and thus exert a vibra 
tion danping influence on the rotor. 

Self-excited vibration phenomena are characterized by an available 
st^jply of energy and by a zero or negatively danped system. Here 


System = Mechanical structure plus Working fluid 

Energy source = Mechanical sources plus Working fluid 

Ihe rotor/bearing/pedestal subsystem, including the hydrodynamic bearing 
film, is generally positively danped. But the ability of bearings to dissi- 
pate energy is a function both of their design (style, dimensions, operat- 
ing conditions, etc.) and their location along the rotor axis. Thus, posi- 
tive damping of the subsystem cannot be assured a priori. Including the 
working fluid will either further increase or decrease the system damping. 
If overall damping capacity exceeds excitation, energy will flow out of “me 
system and stability of the rotor will follow or be assured. But if exci- 
tation exceeds the capacity of the system to dissipate the energy, the net 
influx or accumulation of energy will rapidly produce rotor instability. 


It is probably desirable to reduce the con^pt of energy flux with 
respect to each corrponent of the system to a danping c^fficient 6 appro- 
priate to each conponent. The system stability criterion may then be 
expressed as 


total '^fluid '^mech ^ 


( 1 ) 


Use of an energy criterion for evaluating turboroachinery system stability 
was suggested by Carta (1967, ref. 1) and further explained by him (ref. 2) . 
A description of its application in a design procedure for axial flow com- 
pressor blading was given by Mikolajczak (1975, ref. 3) and a similar ap- 
plication in England was outlined by Halliwell (1977, ref. 4). In 1978, 
Tompson (ref. 5) introduced the stability incronent for multistage centrif- 
ugal cotpressors which, vpon ccmparison, proves to be identical to the cri 
terion of eq. (1) . 


NATURE OF EXCITATION 

The unbalanced force which arises fron the perturbed pressure and 
which is applied to the rotor and leads either to excitation or damping of 
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its vibration/ can be resolved into direct and cross— coijpled oonponents. 

For rotor shaft displaaement dependent forces, such as developed in shaft 
end or blade shroud labyrinth seals, the resolution is conveniently taken 
along and perpendicular, respectively, to the line between the bearing and 
journal centers. For rotor shaft orbital veZoaity dependent forces, such 
^ developed by the working fluid in a turbomachine stage, the resolution 
is taken along and perpendicular, respectively, to a virtual radius to 
which the orbital velocity is perpendicular. The vector diagram of these 
forces is shown in fig. 1. 

In each case, the conponent of the unbalanced force acting along the 
line of centers or along the virtual radius is denoted the direct ooipon— 
snt. If it is in the direction of n 2 or n 2 *, it has small stiffening 
effect on the system, somewhat raising the critical speeds and iirproving 
the st^ility. If the direct conponent acts in a negative n 2 or n 2 * 
direction, it has a sli^t softening effect, reducing critical speeds and 
stabilit^?^ margin. 

In conparison, the conponents of the unbalanced force acting perpen- 
dicular to the line of centers or to the virtual radius are the cross- 
coip>led ccitponents and directly affect the flux of energy to or from the 
rotor. If the cross- coupled conponent is in the ni or ni ' direction, it 
adds energy, excites the rotor, and is a destabilizing influence. If it 
acts in a negative m or m ' direction, is ranoves energy from the rotor 
and is a stabilizing influence. The conception of the effect on rotor 
system stability of the cross-coupled force catponents was first intro- 
duced by Kapitza (1939, ref. 6). 

The unbalanced force, which affects both rotor vibration anplitude 
and system stability, has been attributed to perturbations of the pressure 
on the turbomachine stage passage surfaces. In particular in a reaction 
turbine stage, the effect in the rotor blade passage is conprised of the 
effect of the rotor orbital velocity and the influence of the seal leak- 
age flew in modifying the throughflcw velocity conditions. In addition 
thsre is the effect of the nonuniform circumferential static pressure due 
to nonuniform seal leakage around the periphery of the rotor blade shroud. 
Both passage and seal leakage effects are induced by rotor vibratory motion. 
The conception of the effect on rotorctynamic stability of the seal leakage 
^ it affects both blade performance cind generates an unbalanced force was 
first introduced by Thomas (1958, ref. 7) 

In more recent times, some valuable expositions have been made of the 
concepts given above. Ehrich (1972, ref. 8; 1973, ref. 9) identified an 
aeroelastic tip-clearance effect and presented the associated destabilizing 
force vector diagram. Alford (1965, ref. 10) in an oft-cited paper, iden- 
tified two aerod^amic disturbing forces, the circumferential pressure 
variation in the seals and the variation of local efficiency in the rotor 
blade flow processes, where both were attributed to blade tip clearance. 
Based on very astute interpretation of the ejqperimsntal test results of 
four gas turbine engines, Alford also concluded that rotor whirl resulting 
fretn fluid c^mamic excitation was in the direction of rotor rotation, that 
vbirl anplitude increased both with increased power output (i.e. , in large 
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part with an increase in mass flow rate) and with decreased inlet tertper- 
ature (i.e., with increased inlet density). These observations seem to be 
confirmed by ejperience accumulated up to the present time. 

Pollman, Schwerdtfeger, and Termuehlen (1978, ref. 11) provided an 
extensive account of excitation mechanisms investigated by a European steam 
turbine manufacturer. Excitation of a fluid d^amic nature was ccnprised 
of the effects of seal leakage flow and load dependent influences, so-called 
"steam whirl." The latter were clearly defined as modification of local 
rotor blade passage efficiency caused by steam leakage flow and resulting 
in a local variation in torque developed in the blading. An unbalanced 
force in the direction of rotor rotation follcws from the torque variation. 

In Russia, work in this area is carried out at two research institutes, 
the Mdscow Energetics Institute (MEI) and the Central Boiler-Turbine Insti- 
tute (TsKTI) in support of steam turbine manufacturing. Kostyuk, et al. 

(1974, ref. 12) identified three nonoonservative forces leading to self- 
excited vibrations. These were oil film bearing forces, seal leakage 
forces, and circumferentially varying blade forces due to seal leakage. 
Further Kostyuk, et al. (1975, ref. 13) reported work on two additional and 
important effects, that of the influence on the blading of passing through 
the wakes of upstream blades and that of unsteady separation on the blade 
surfaces due to condensation shock wave/boundary layer interaction which 
occurs at transonic and mildly supersonic velocities when the steam condi- 
tion is sijpersaturated. These latter two effects clearly excite blade vi- 
bration, vdiile the former induce rotor vibration. In a later report Kostyuk, 
et al. (1978, ref. 14) presented enperimental results of labyrinth seal force 
measurements. Olinpiev (1978, ref. 15) presented alternative work on vary- 
ing blade forces, on shroud band seal leakage forces, and on labyrinth seal 
forces. He also mentioned that seal configurations could be devised to 
convert the influence from excitation to damping and that the circumfer- 
entially varying blade forces taken prcportional to the rotor displacoient 
producing them lead to an effective rotor stiffness coefficient. 

NECESSITY OF UNSTEADY FLOW 

Dean (1959, ref. 16) has explained the paradox encountered in applying 
E\jler's equation (the equation of motion) to the problem of describing the 
change cf state of the working fluid as work is transferred in a turbo— ^ 
irHchine. Consider a frictionless, reversible, adiabatic process in an ideal- 
ized machine. This is simply an inviscid fluid subject to an isentropic 
process, which assumptions do not, in the present argrment, invalidate the 
conclusions. The energy equation is simply and correctly given as 

-Wg = ho 2 “ ho 1 

For a compressor, ho 2 ^ ho i but for a turbine ho 2 ho 1 and the sign of the 
result is entirely consistent with the conventions of thermodynamics. 

However, Euler's equation (the differential form of the equations of 
motion) and Bernoulli's equation (integral form of the eqxrations of motion) 
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are not applicable in describing the process in detail whereby the shaft 
work changes the state of the working fluid as it passes throu^ blading of 
turbanachines . Instead, as Dean shews in detail for compressible flow, the 
equations of motion become 

Eho/Dt = (1/p) 9p/9t (3 

and for inooitpressible flow 

Dpo/Dt = 9p/9t (4 

where the material derivative is the operator having unsteady and convective 
parts: 

D/Dt = 9/9 t + (c.V) 

■^us, while the energy equation (2) dotonstrates that stagnation enthalpy 
(or stagnation pressure) is changed as a consequence of shaft work, eq. (3) 
or eq. (4) shows, first, that the changes are with respect both bo time and 
to changing location in the fluid field and, second, that both are a func- 
tion of the time- varying or unsteady pressure fluctuations. Therefore, 
vork transfer in a turbcmachine requires am unsteat^ flow process. 

Spannhake (1930, ref- 17) perhaps was the first to point this out to turbo- 
machine designers and analysts. 


Three excellent reviews of unsteady flow phenomena can be identified, 
v4iere in turn a large body of literature is cited. Mikolajezak (1976, ref. 
18) pointed out that unsteaefy effects affect aerodynamic performance, 
aeroelastic and rotorc^amic performance, and the generation of noise. 
Platzer (1978, ref. 19) provided an exhaustive account of unsteac^ phencm- 
turbcmachinery including the unsteady effects in machines e>^x 3 sed to 
uniform flew, the response to distorted inlet flew, surge and rotating stall 
conditions, as well as blade flutter. He categorized various prediction 
methods for these phenomena and described ej^rimental studies. McCroskey 
(1977, ref. 20) included unsteac^ turbomachinery flew effects in a general 
review of unsteaefy fluid cfynamics. Even so, he was able to make critical 
comments about flutter, inlet distortion, unsteac^ transonic flow, rotating 
stall, transito:^ diffusor stall, surge, vortex shedding frem bluff bodies, 
^d the def ini toon of an unsteac^ Kutta-Joukewski condition. Finally, he 
included a projection of research needs and future developments. 

Two tome scales occur in the description of unsteac^ flow phenomena, 
the particle transport time: 

Characteristic length/Fluid velocity = I/w [tine] 

and the period of the pressure fluctuation: 

l/Fluctuatoon frequency = 1/v [time] 

Their ratio becemes 

(Lv)A? = k, the reduced frequency ( 5 ) 

When the particle transport time is short with respect to the fluctuation 
period, k ->■ 0 and the flew is said to be quasi-stea(^. Practically, for 
k 'X/ 0.1, the approximation of stea<^-flow analysis is often acceptable. 
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Lastly, unsteady phenomena iray be transient or periodic. While a 
single transient event my have catastrophic consequences for a tiurbo- 
itachine, such events occur over a relatively long period of time and are 
usually under the control of the process control system. On the other 
hand, periodic events characterize self— excited vibrations and are of 
concern here. Such periodic events fall into three categories ^ a f motion 
of both their frequency and a characteristic length of the p>articipating 
fluid enserrble. For a characteristic length of the order of the blade 
thickness, the blade passing frequency is appropriate. Frequencies of the 
order of the shaft rotational speed are associated with phenomena of the 
size of the stage passages, while phenomena of the size of the bearing span 
cxccur at frequencies of the order of the overall duct system resonance. 
Vibration exciting mechanisms in each of the categories will be described 
in the next section. 

SIZE/FREQUENCY CLASSIFICATION OF PHENOMENA 

Three classifications according to size and frequency of the msteady 
flow phenonena have been noted above. In each classification one finds a 
nuitber of phenomena each one of v^iich inherently possesses a mechanism for 
exciting vibration of the associated structural elonent. At the order of 
the blade passing frequency there occur axial compressor and axial turbine 
blade flutter together with blade vibration due to cutting upstream blade 
wakes as well as condensation shock wave oscillations. The last condition 
occurs in steam turbines v\iiere the ej^ansion crosses the saturation line 
and the flow is at transonic speed in order to acccmtnodate the large steam 
volume within a reasonable passage size. The frequency of the shock wave 
oscillations is below the blade passing frequency but much above the shaft 
rotational speed so it is included in this category. 

At the order of the shaft rotational speed, inlet flew distortion, 
gusts in the flow and rotating stall in the rotor blading or the stationary 
interstage passages are external influences ntposed on the rotor. In 
addition rotor vibration due to, say, residml mass unbalance is capable 
of inducing both displacerient dependent and orbital velocity dependent ex 
citation of the rotor. Finally at the order of the system resonance fre- 
quency, surge can occur. 

All of these topics are listed in table 1 together with a group of 
references for each phenomena. The references have been chosen frem among 
a very large number because each in turn identifies and to a great extent 
evaluates the contributions of a subsequent set of references. In comec- 
tion with table 1, the writer is indebted to William G. Steltz, Fewer 
GeneraUon Divisions, WOstinghouse Electric Corp. for the references on 
condensation shock wave phencmena. 

Consistent with the present workshop on rotordynamic instability, 
phenomena at the order of the shaft synchronous frequency will be enphasized 
and blade vibration and system surge will not be considered further. Fur 
ther we mderstand that from the msteady flew phenomena there needs to be 
derived the amplitude, frequency and phase angle of a forcing fmetion 
vbich, for analytical purposes, can then be imposed on the rotor system. 
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But it cannot be expected that all topics in table 1 are so well developed 
fluid dynamically that clear, realisUc and practical forcing functions 
will be available for our use. 

ORBITAL VELOICITY-DEPENDEISIT EXCITATION 


The orbital velocity induces the unbalanced fluid dynamic force v^ich 
IS imposed on the rotor of a turbcxnachine. In both figures 2 and 3, the 
rotor centerline describes a trajectory or orbit during lateral vibratory 
motion. Uie instantaneous tangential velocity of the shaft center along 
the trajectory is the rotor orbital velocity. 

The angular velocity of the rotor results in a tangential velocity at 
each location on the rotor proportional to its radius frcxn the shaft center. 
Then for toe leading edges of the blades as shown in figures 2 and 3, \toen 
toe velocity due to rotation and the orbital velocity are in the same sense, 
they add and in opposite sense, they subtract. The effect is as if 

the rotor rotation took place instantaneously about a virtual center dis- 
placed from the geometric shaft center as shown in the figures. Those flow 
passages at the greater radii from the virtual center handle a somewhat 
greater mass flew and exchange a somewhat greater angular momentum with the 
fltod toan do those at the smaller radii. Blade surface pressure dis- 
tributions in toe flow passages will now vary with the virtual radii, not 
the geometric radii, and will be nonuniform from passage to passage. Sum- 
mation of toe pressure distributions in all passages at an axial location 
ori the rotor will thus result in an unbalanced force on the rotor. 

In the case of toe radial conpressor in fig. 2, the rotor orbital 
motion causes the blade leading edges to move in and out radially and, as 
the orbital velocity first adds then subtracts fresnn the rotational velocity, 
to accelerate then decelerate tangentially with respect to a steady fluid 
flow into the passages. Thus, relative to the rotor, the blade appears to 
be moving through a gust having both lateral and longitudinal components of 
perturbed velocity and pressure. Such perturbations move in a series of 
waves over the length of the blade with a frequency equal to the rotor 
lateral vibratory frequency. In practice this frequency is the sub- 
synchronous frequency of the rotor/bearing/pedestal systan. Isay (1958, 
ref. 69) appears to offer toe only known analysis of such unsteady flow. 

to much the same way, orbital motion induces velocity perturbations, 
mcxlified incidence angles, etc. in an axial conpressor rotor as shown in 
fiff* 3. to addition the same orbital motion produces varying blade tip 
clearance or blade shroud band clearance about the periphery at a given 
rotor axial lotion yielding varying leakage. The local variation of the 
leakage modifies the local pressure rise in an axial cenpressor stage and 
modifies toe local pressure drop in a turbine stage. Thus the velocity 
perturbation frean the orbital motion and the pressure pjerturbation from 
the rotor orbital displacanent combine to perturb the passage flow non- 
uniformly fron passage to passage resulting in an unbalanced force on the 

In the case of the axial flew steam turbine, the cembined influence 
of leakage and orbital notion is called steam vhirl. 
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EXAMPLE: THE CENTRIFUGAL CCMPRESSOR STAGE 

At the outset of this work, ref- 69 was unknown to the writer and 
circumstances did not permit developing the unsteady flew analysis. In- 
voking quasi-steady conditions, a "snap shot" was made enccxipassing the 
flow in every inpeller channel where each is perturbed differently by the 
shaft orbit velocity. When the blade or channel length is a sufficiently 
small portion of the disturbance wave length, the flow conditions inposed 
for an instant on the blade can be considered approximately steady. 
pressed in terms of the reduced frequency introduced above, k < 0.2 will 
usually be adequate for the approximation to hold. The fluctuation fre- 
quency is given by the shaft subsynchronous frequency. 

A conputer program was organized to analyze each and every impeller 
passage under perturbed, quasi-steady conditions; to resolve the summation 
of the pressure distribution into cottponents of Fi and F 2 contributed by 
each passage; and finally to sum the contributions to Fi and F 2 . Practical 
considerations showed that typical orbital velocities are 0.05 to 0.20 m/s 
(2 to 8 in./s) in conparison bo the passage relative fluid velocity of 46 
to 185 it/s (150 to 600 ft/sec). The unbalanced force resulting from such 
a small perturbation of the velocity proved, upon calculation, to be linear 
in the orbital velocity. Thus F^ = b^ F 2 = b 2 u^, and the contribution 
of the working fluid to the system stability criterion becanes 

^fluid = 

The negative sign is needed because a positive excitation corresponds to a 
negative damping. For an extended explanation of these ideas as well as 
numerical results for a cempressor ccxrprising fourteen stages in two casings, 
the reader is referred to ref. 5. 

CONCLUDING REMARKS 

The system stability criterion given in ref. 3 is shewn in eq. 1. For 
the relatively siitple blade structural system considered in ref. 3, use of 
the ideal, zero boundary between stability and instability is evidently 
justified. Such is not the case in the much more oonplex rotor/bearing/ 
pedestal/working fluid system. As a result of analyzing 115 individual 
centrifugal coitpressor stages v^uLch were contained in 20 multistage bodies, 
the average coefficient per stage needed to assure stability may be suggest- 
ed as 1-85. Referring again to eq. 1, only is the si±)ject 

of this paper \diile must be determined by other rotodynamic analysis. 

Each length/frequency classification of excitation phenomena, once 
analyzed, will admit suggestions for designing in additional damping of 
a fluid dynamde nature. Regarding blade flutter, ref. 3 suggests that an 
entire blade fluid and structural design system has been organized and can 
yield optimized results. Regarding distorted inlet flew to axial compress- 
ors, criteria are emerging (refs. 70, 71) \diich will lead to configurations 
vhich are insensitive to such external influences. Work in this direction 
for radial turbomachines lags that of axial machines. Finally overall 
system design ignores completely the opportunity to add damping to minimize 
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the effects of surge. 


Regrettably the limitation on scope of this paper due to both prepar- 
ation time and space has precluded including examples of other phenon^a 
in table 1, an evaluation of the numerical fluid mechanical methods v^ich 
are used to obtain quantitative predictions of the damping coefficients 
and evaluation of various experimental techniques for observing unsteady 
flow and deriving the force and damping coefficients from such observations. 

SYMBOLS 


ai , a.2 


ai f a.2 

bl/ b2 


c 


F 

1 



F 

2 



h 

ho 

k 

L 


ni, n2/ n3 


P 

Po 

r 

u 

u. 


w 


displacement dependent fluid force coefficient, N/m 
(Ibf/in. ) 

stator-fixed coordinates syston, m (ft) 

orbital velocity dependent flioid force coefficient, 
t^s/m (Ibf-s/in.) 

fluid absolute velocity, rti/s (ft/s) 

cross-coupled force ccxrponent due to fluid dynamic 
excitation, N (Ihf ) 

direct force ccxrponent due to fluid dynamic excitation, 
N(lbf) 

static enthalpy, m-NAg(ft-Ibf/slug) 
stagnation enthalpy, m-NAg (ft-lbf/slug) 
reduced frequency 
characteristic length, m(ft) 
rotor-fixed coordinate system, m(ft) 
static pressure, kPa(psia) 
stagnation pressure, kPa(psia) 
radial ccordinate, m(ft) 
blade peripheral velocity, m/s (ft/s) 
rotor orbital velocity, m/s (in./s) 
fluid relative velocity, m/s (ft/s) 
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shaft work, m-NAg (f t-lbf/slug) 
damping coefficient, N-s/m (Ibf-s/in. ) 
rotor displaconent, m(in.) 

fluctuation frequency (rotor subsynchronous frequency) , 

1/s 

density, kg/m^ {slugs/ft ) 
angular velocity, rad/s 
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Table 1 Size-frequency Classification of Vibration Excitation Mechanisms 


Topics 


References 


Blade thickness/blade passing phenomena 





Axial conpressor blade flutter 

19, 

3, 2, 21 


Axial turbine blade flutter 

22 




Wake flow excitation 

23, 

24, 25, 

26, 

27 

Condensation shock wave excitation 

28, 

29 , 30 , 

31, 

32 

Passage dimens ion/sha ft synchronous phencxnena 





Axial configuration: 





Rotor displacement dependent 

7, 

10, 33, 

34 


Orbital velocity dependent 





Inlet distortion 

35, 

36, 37, 

38, 

39, 


40 




Gusts 

41, 

42, 43, 

44, 

45 

Rotating stall 

! 46, 

47, 48 



Radial configuration: 





Rotor displacement dependent 

49, 

50, 51, 

52 


*Orbital velocity dependent 

; 5 




Rotating stall 

’ 53, 

54, 55, 

56 


Bearing span/systems resonance 





Sixrge: General 

57, 

58 



Punp 

59, 

60, 61 



Conpressor 

62, 

63, 64, 

65, 

66, 


67, 

68 




* An example will be presented in tbe text. 
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F = F + F, 
a b 


^2/ 


T\2 


direct cotponent unit vectors 


Where F is the displaconent dependent force, 

3 . 

F = F n - F n 
a ai 1 ^2 2 

and Fj^ is the velocity dependent force. 


Fk = F, n' 
b bi 1 


F, n^ 

02 2 


Further the rotor displacement vector 6 


6 = 6{-n^) 


and the orbital velocity vector is 

since F^ and Fj^ are linear in 6 and respectively 

^bl = ''b2 = 

(Note: F^= Fj^^, F^ = text) 


Figure 1 Description of Fluid Dynamic Forces 
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Figure 2 - Fluid Dynamic Excitation Force of a Radial Conpressor Impeller 
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Figure 3 Fluid Dynamic Excitation Force of an Axial Conpressor Rotor 






